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ABSTRACT 


Experimental  investigations  were  conducted  in  a  Mach  3 
combustion  tunnel  to  determine  heat  release  characteristics  of 
a  hydrogen-air  mixture  after  it  had  passed  through  a  normal 
shock.  Inlet  temperature  was  varied  from  1800  to  3000°R. 
Heat  release  was  indicated  by  three  types  of  measurements: 

1.  Combustion  and  total  temperature  rise  calcu¬ 
lated  from  gas  analysis 

2.  Total  pressure  loss  caused  by  heat  addition 

3.  Ultra-violet  emission  at  the  OH  emission 
frequencies,  indicating  an  H2-O2  reaction 

Combustion  efficiency  was  found  to  increase  with  temper¬ 
ature  and  to  be  independent  of  fuel- air  ratio. 

Fuel  passing  through  an  oblique  shock  was  observed  to 
emit  a  radiation  similar  to  that  observed  from  the  normal 
shock  combustion. 
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NOMENCLATURE 


at  Speed  of  sound  at  total  temperature,  ft/ sec 

^2  Excess  fuel  over  stoichiometric  equivalence  ratio 

E  Global  activation  energy,  cal/gm  mole 

g  Dimensional  constant,  Ibm  ft/lbf  sec^ 

h  Enthalpy,  Btu/lbm  mole 

M  Mach  number 

m  Molecular  weight 

mi  Mole  fraction 

p  Pressure,  psia 

Prt  Turbulent  Prandtl  number 

Pt  Total  pressure  behind  a  normal  shock,  psia 

Pto  Plenum  pressure,  psia 

Qj  Heat  of  reaction,  Btu/lbm  mole 

Ro  Universal  gas  constant,  energy/lbm  mole  °R 

T  Temperature,  °R 

Tto  Total  temperature  with  e  =  0,  “R 

V  Velocity,  ft /sec 

X  Axial  distance  behind  shock  wave,  in. 

Y  Vertical  distance  from  tunnel  centerline,  in. 

Ti  Preheater  fuel  concentration  in  equivalence  ratio 

Test  fuel  concentration  in  equivalence  ratio 
Z  Horizontal  distance  from  tunnel  centerline,  in. 

y  Specific  heat  ratio 

f  Combustion  efficiency  of  test  fuel 

T/  Overall  combustion  efficiency,  test  fuel  and  preheater  tuei 

fj  Density,  Ib^/ft^ 

r  Time,  microsec 

[A]  Moles  of  component  A  per  unit  volume 
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(A)  Volumetric  fraction  of  component  A  in  mixture  (dry) 

SUBSCRIPTS 

c  Centerline 

f  Test  fuel  only 

p  Preheater  only 

pf  Preheater  and  test  fuel 

t  Total 

1  Initial  conditions 

2  Final  conditions 
Free  stream 


xii 


AEDC-TDR. 62-78 


1.0  INTRODUCTION 


Although  the  chemical  reaction  rate  of  hydrogen  has  been  studied 
for  many  years,  little  information  is  available  on  the  effect  of  the  aero- 
thermodynamics  of  a  supersonic  airstream  on  the  chemical  reaction 
rate.  Many  theoretical  analyses  of  devices  using  supersonic  combustion 
have  been  made  -  for  example  ramjet  engines  (Refs.  1,  2,  and  3).  How¬ 
ever,  the  problem  of  how  this  combustion  can  be  stabilized  and  controlled 
has  not  been  solved. 

Although  detonations  in  shock  tubes  have  been  studied  previously, 
the  attempt  to  work  with  "standing"  waves  is  relatively  recent.  Nicholls 
(Refs.  4  and  5)  used  a  normal  shock  wave  produced  by  a  highly  under¬ 
expanded  free  jet  nozzle.  Gross  (Ref.  6)  and  Rhodes  and  Chriss  (Ref.  7) 
used  a  normal  shock  produced  by  the  intersection  of  two  wedge  shocks. 

This  report  discusses  the  work  done  in  the  Rocket  Test  Facility  (RTF), 
Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Systems 
Command  (AFSC),  using  the  Supersonic  Combustion  Tunnel  (SCT)  as  a 
continuation  of  the  work  reported  in  Ref.  7.  Extensive  data  are  pre¬ 
sented  on  the  flow  field  characteristics  in  the  zone  of  the  two- wedge  shock 
and  on  the  chemical  reaction  rate  of  hydrogen- air  mixtures  in  this  flow 
field  when  the  temperature  level  is  high  enough  for  chemical  action  to 
occur.  The  work  was  directed  toward  establishing  a  simple  aerodynamic 
model  in  which  combustion  phenomena  could  be  studied  with  uncompli¬ 
cated  instrumentation  which  could  measure  conditions  at  a  point  in  the 
stream.  In  these  experiments,  the  objective  was  to  verify  whether  heat 
release  actually  occurred,  and,  if  so,  to  what  extent.  This  was  done  by; 

1.  Sampling  combustion  gases  through  a  probe  specially  designed  to 
promote  rapid  quench  of  the  chemical  reactions  and  determining 
total  temperature  and  other  properties  of  the  gas  from  an  analysis 
of  the  gas  composition; 

2.  Making  extensive  pressure  surveys  and  comparing  the  measure¬ 
ments  to  those  predicted  by  theoretical  equations  to  determine 
whether  deviations  were  caused  by  heat  release;  and 

3.  Measuring  the  radiant  emission  from  the  reaction  with  a  spectrom¬ 
eter  to  determine  whether  radiation  was  caused  by  OH  ion  emission, 
and  hence  that  an  H2  -  O2  reaction  was  in  progress. 

Ultimately,  knowledge  gathered  from  these  studies  could  be  useful  in 
the  design  of  high-speed  propulsion  systems. 


Manuscript  released  by  authors  March  1962. 
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2.0  APPARATUS  AND  PROCEDURE 


2.1  SUPERSONIC  COMBUSTION  TUNNEL 

A  detailed  description  of  the  Supersonic  Combustion  Tunnel  (Fig.  1) 
is  presented  in  Ref.  7.  Air,  preheated  by  an  indirect-fired  heater  to 
1460“R,  enters  the  plenum  where  the  core  of  the  nozzle  airflow  is  further 
heated  to  a  maximum  of  SSOO^R  by  the  combustion  of  hydrogen.  The 
nozzle  is  two-dimensional  and  initially  had  an  exit  area  3  in.  wide  by 
5  in.  high.  It  is  operated  as  a  free  jet  with  a  discharge  Mach  number  of 
3.  1  inside  the  test  rhombus. 

During  these  tests  the  tunnel  width  was  increased  to  six  inches.  This 
did  not  seem  to  have  any  effect  on  the  flow  parameters  except  to  widen 
the  zone  in  the  center  which  was  free  of  disturbances  from  the  sidewalls. 

The  procedure  for  establishing  the  desired  conditions  in  the  test 
section  is  described  in  detail  in  Ref.  7. 


2.2  TEST  FUEL  INJECTION 

A  decision  was  made  early  in  the  program  to  use  a  thin  wedge  fuel 
injector  installed  in  the  supersonic  zone  of  the  nozzle  (Figs.  1  and  2), 
rather  than  the  tube  fuel  injector  near  the  tunnel  throat  used  in  the  tests 
reported  in  Ref.  7.  The  purpose  of  this  modification  was  to  reduce  the 
possibility  of  partial  burning  of  fuel  near  the  injector  where  the  static 
temperature  was  quite  high.  Some  experiments  with  throat  fuel  injection 
at  the  higher  temperatures  indicated  that  this  type  of  partial  burning  was 
occurring.  A  photograph  of  the  fuel  injector  is  shown  in  Fig.  3.  Further 
analysis  of  the  wedge  aerodynamics  will  be  found  in  section  3.  1.  1. 


2.3  WEDGES 

Two  types  of  wedges  were  used  to  produce  shock  systems.  The  first 
was  a  single,  water-cooled,  30-deg  wedge  for  producing  an  oblique  shock 
wave.  This  wedge  was  mounted  in  a  window  frame  as  shown  in  Fig.  4. 
The  second  type  consisted  of  a  set  of  two  25-deg  wedges  used  to  produce  a 
pair  of  oblique  shocks  and  a  connecting  normal  shock.  The  shock  system 
produced  and  the  location  of  the  wedges  are  shown  in  Fig.  2.  The  photo¬ 
graph  in  Fig.  1  shows  the  wedges  retracted  from  the  test  position. 
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2.4  INSTRUMENTATION 

All  temperatures  except  the  test  section  total  temperature  were 
measured  with  chromel-alumel  thermocouples  and  recorded  on  a  record¬ 
ing  potentiometer.  The  total  temperature  in  the  test  section  was  cal¬ 
culated  using  the  analysis  of  the  gas  as  shown  in  Appendix  A. 

All  air  pressures  except  the  total  and  static  pressures  in  the  test 
section  were  measured  with  a  120-in. ,  mercury-filled,  manometer  board 
and  recorded  photographically. 

Hydrogen  pressures  were  measured  with  transducers  and  recorded 
on  a  recording  potentiometer. 

A  water-cooled  probe  was  used  alternately  for  obtaining  gas  samples 
and  total  pressure  measurements  in  the  test  section.  When  required,  a 
second  probe  was  used  to  measure  the  static  pressure.  Figure  5  shows 
probes  with  and  without  the  static  tube.  The  total  pressure  and  static 
pressure  were  recorded  as  a  function  of  position  (the  Z  tunnel  coordinate) 
on  an  x-y  plotter.  A  schematic  diagram  of  the  total  pressure,  static 
pressure,  and  gas  sampling  system  is  shown  in  Fig.  6. 

Samples  were  pumped  from  the  probe  with  a  dry  diaphragm  pump 
which  maintained  a  pressure  of  2  to  3  psi  inside  the  probe  and  delivered 
the  samples  to  the  analyzers  at  slightly  above  atmospheric  pressure.  A 
water  separator  and  chemical  drier  were  used  in  the  sampling  lines  and 
all  analyses  were  made  on  a  dry  basis.  Local  concentration  of  hydrogen 
in  the  test  section  was  measured  by  a  thermal  conductivity  meter  de¬ 
scribed  in  Ref.  7.  Local  concentration  of  oxygen  was  measured  by  a  meter 
which  detects  the  change  in  the  magnetic  susceptibility  of  the  gas.  The 
oxygen  meter  was  calibrated  using  nitrogen  for  a  zero  and  atmospheric 
oxygen  for  20.  95-percent  oxygen.  The  oxygen  and  hydrogen  concentrations 
were  recorded  on  a  recording  potentiometer.  An  analysis  of  the  mecha¬ 
nism  of  sampling  is  discussed  in  Appendix  B. 

Shock  position,  axial  and  vertical  position  of  the  test  section  probe, 
and  location  of  emission  from  high  temperature  zones  were  determined 
photographically.  The  test  section  was  photographed  with  a  schlieren 
camera  using  either  collimated  light  from  a  spark  source  or  direct  emis¬ 
sion  from  the  test  section.  The  schlieren  spark  was  collimated  to  a 
6- inch  beam  by  an  aerial  camera  lens.  This  light  and  the  emission  from 
the  test  section  were  focused  on  a  4  in.  by  5  in.  camera. 
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2.5  DATA  REDUCTION 

Pressure  data  were  read  from  manometer  board  photographs  or 
X  -  y  recorder  charts  and  reduced  manually.  Temperature  and  composi¬ 
tions  were  obtained  from  a  computer  solution  of  the  equations  given  in 
Appendix  A.  An  estimate  of  the  accuracy  of  the  raw  and  reduced  data  is 
also  given  in  Appendix  A. 


3.0  RESULTS  AND  DISCUSSION 


Shock- induced  combustion  occurs  when  a  premixed  combustible 
mixture  is  heated  to  above  its  ignition  temperature  by  passing  through  a 
shock  wave.  The  most  familiar  example  of  this  is  the  classical  detonation 
wave  in  which  a  constant  area  boundary  condition  is  imposed  on  the  flow. 

In  this  case  the  energy  needed  to  drive  the  shock  wave  is  supplied  by  the 
combustion.  This  is  the  type  of  phenomenon  seen  in  a  detonation  tube. 

Here,  the  final  conditions  are  determined  by  the  initial  conditions  and  the 
application  of  conservation  of  mass,  momentum,  and  energy. 

Previous  workers  with  shock- induced  combustion  waves  have  assumed 
that  the  thermodynamics  of  constant  area  classical  detonations  are  appli¬ 
cable  to  all  shock- induced  combustion  waves  (Refs.  4  and  7).  It  would 
seem,  however,  that  other  aerodynamic  processes  for  releasing  heat  using 
shock- induced  combustion  would  be  possible  if  different  boundary  conditions 
are  applied.  It  is  possible  to  describe  theoretically  a  constant  pressure, 
shock- induced  combustion  wave  where  the  gases  heated  by  the  combustion 
expand  so  that  the  static  pressure  in  the  combustion  zone  is  the  same  as  the 
static  pressure  behind  the  shock  wave  without  combustion.  Physically 
this  might  occur  if  a  combustible  core  of  gas  passing  through  a  shock  wave 
had  a  relatively  small  diameter  compared  to  the  reaction  zone  length  so  that 
the  expansion  of  the  heated  gases  would  not  appreciably  change  the  static 
pressure  behind  the  shock  wave.  It  might  also  occur  in  a  variable  area  duct 
where  the  contour  exactly  matched  the  rate  of  heat  addition  to  maintain  con¬ 
stant  static  pressure.  As  long  as  the  rate  of  expansion  caused  by  heat 
addition  is  low  compared  to  the  velocity  of  the  burning  gas,  there  will  be 
negligible  radial  pressure  gradients  across  the  combustion  zone;  and,  if 
a  step  profile  can  be  assumed,  a  one -dimensional  flow  approximation  will 
describe  the  flow.  As  the  ratio  of  the  diameter  to  the  length  of  the  reaction 
zone  increases,  the  pressure  at  the  center  of  the  reaction  zone  approaches 
the  constant  area  value  because  three-dimensional  expansion  will  cause 
sizeable  pressure  gradients  from  the  center  to  the  outside  of  the  combustion 
zone.  Constant  pressure  shock- induced  combustion  is  then  limited  to  com¬ 
bustion  zones  with  long  enough  reaction  times  and  slow  enough  expansion  so 
that  these  three-dimensional  effects  are  negligible.  This  type  of  reaction 
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cannot  produce  a  normal  shock  wave  ahead  of  it  because  the  same  size 
and  reaction  rate  criteria  which  allowed  the  one- dimensional  approxima¬ 
tion  will  preclude  the  pressure  rise  necessary  to  produce  the  shock  wave. 
Again  the  terminal  conditions  of  the  constant  pressure  case  can  be  deter¬ 
mined  explicitly  from  the  inlet  conditions  and  the  continuity  equations. 

A  more  complicated  case  arises  when  there  is  an  axial  static  pres¬ 
sure  gradient  behind  the  shock  wave.  For  a  small  diameter  combustion 
zone,  this  pressure  can  exist  independent  of  the  combustion.  If  the  com¬ 
bustion  zone  size  and  the  reaction  rate  are  small  enough  and  step  profiles 
are  assumed,  a  one -dimensional  approximation  will  in  this  case  ade¬ 
quately  define  the  flow;  however,  the  final  conditions  will  depend  on  the 
velocity  at  which  the  heat  is  added  as  well  as  on  the  inlet  conditions.  The 
terminal  conditions  are  not  defined  by  the  inlet  conditions  and  the  total 
heat  release  as  they  are  in  the  other  cases  mentioned.  In  this  case  the 
static  temperature  and  therefore  the  rate  of  the  reaction  will  be  influenced 
by  the  pressure  gradient. 

The  last  case  is  the  model  which  was  used  for  the  analysis  of  the  data 
from  the  SCT.  A  static  pressure  gradient  exists  behind  the  shock  system 
which  depends  on  the  interaction  of  the  shock  and  expansion  waves.  The 
assumption  that  this  pressure  gradient  is  independent  of  combustion  was 
verified  by  the  data.  Total  pressure  loss  for  a  given  amount  of  heat  re¬ 
lease  depends  on  the  location  in  this  pressure  field  where  the  heat  is  re¬ 
leased. 

All  the  possible  boundary  conditions  are  not  limited  to  those  mentioned 
thus  far,  since  by  application  of  the  proper  pressure  gradient  or  area 
change  it  would  be  theoretically  possible  to  achieve  heat  release  at  con¬ 
stant  Mach  number  or  constant  static  temperature,  although  these  might 
be  impossible  to  achieve  in  practice.  In  an  actual  combustion  system, 
combinations  of  dtffusional  and  shock-induced  combustion  may  exist,  and 
one- dimensional  theory  would  not  be  adequate.  This  system  would  be  even 
more  difficult  to  analyze,  although  over-all  performance  may  be  measurable. 


3.1  FLOW  FIELD  DESCRIPTION 

I 

In  this  section  are  described  the  aerodynamic  and  thermal  conditions 
which  exist  in  the  test  section  of  the  SCT  which  establish  the  boundary 
conditions  for  the  combustion.  This  information  is  based  on  total  pres¬ 
sure  and  total  temperature  traverses  and  on  a  limited  amount  of  static 
pressure  data. 
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A  general  schematic  of  a  Mach  reflected  normal  shock  system  is 
shown  in  Fig.  7.  The  "normal"-shock  wave  is  produced  because  the 
turning  angle  necessary  to  straighten  the  flow  behind  the  first  oblique 
shock  is  greater  than  the  maximum  turning  angle  at  this  Mach  number, 
so  that  a  simple  reflection  cannot  occur.  The  flow  at  the  shock  inter¬ 
section  is  not  turned  parallel  and  contracts  with  supersonic  flow  out¬ 
side  and  subsonic  flow  inside  the  slip  line.  Mixing  occurs  along  this 
boundary  resulting  in  the  reacceleration  of  the  center  flow  both  by  con¬ 
traction  of  the  flow  and  by  the  mixing  with  higher  energy  air. 

The  configuration  used  for  the  normal  shock  studies  has  been  de¬ 
scribed  in  the  section  on  apparatus  and  shown,  with  the  shock  system 
produced  by  this  configuration,  in  Pig.  2. 

3.1.1  Pressure  Effects  Ahead  of  the  Normal  Shock  Wave 


In  the  region  ahead  of  the  normal  shock  wave  the  air  accelerates 
nearly  isentropically  from  the  throat  until  it  passes  around  the  fuel  in¬ 
jection  wedge  which  causes  a  reduction  in  total  pressure.  As  a  result 
of  wedge  interference,  the  total  pressure  and  Mach  number  profiles  of 
the  air  approaching  the  normal  shock  wave  are  distorted  and  have 
minimum  values  at  the  tunnel  horizontal  centerline  (Fig.  7), 

A  calculation  was  made  to  determine  if  the  discrepancy  between  the 
indicated  total  pressures  measured  at  the  centerline  and  the  theoretical 
value  for  a  Mach  3.  1  normal  shock  wave  could  be  explained.  The  pres¬ 
sure  loss  caused  by  the  fuel  injection  wedge  was  compared  to  the  momen¬ 
tum  per  unit  area  at  the  shock  wave  with  the  inlet  momentum  minus  the 
wedge  drag.  It  was  assumed  that  the  dimensionless  velocity  ratio 
(V  -  V^/Vc  -  resulting  from  the  wedge  drag  without  fuel  injection 

is  the  same  as  the  dimensionless  concentration  profile  (y/Xc)  when  there 
was  fuel  injected  (Ref.  8).  A  turbulent  Prandtl  number  of  0.  5  was 
assumed.  Constant  static  pressure  was  assumed  along  the  mixing  zone. 

The  momentum  at  the  shock  wave  was  1/p  /  pV^dV  per  square  inch 
where 

.  y  (V/a,)’  [l  -  (V/a,)‘J'‘  <Ref.  10) 

and  since  for  a  given  test  condition  at  was  constant 

^V/at),t  y  =  (Vc/at  -  yj^t)  hhcKi  y  + 

P  P 

where  Vc/at  was  calculated  from  —  (experimental)  and  -  (isen,  M  =  3.1); 
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V^/at  was  calculated  for  isentropic  M  =  3.1;  and  (7/7^)  was  the  fuel  distri¬ 
bution  at  the  same  total  temperature. 

The  momentum  at  the  fuel  injector  was  assumed  to  be  the  momentum 
calculated  from  isentropic  relationships  at  M  =  3.1  minus  the  wave  drag 
and  skin  drag  on  the  wedge. 

The  following  table  compares  the  momentum  at  the  wedge  with  that 
at  the  shock  for  two  total  temperatures. 


Tt  -  °R 

At  the  Wedge 

pvVp 

At  the  Shock  Wave 

1500 

12.  66 

12.  54 

2750 

12.  11 

12.  17 

The  agreement  between  the  momentum  calculated  at  the  two  stations 
is  close  enough  that  the  reduction  in  total  pressure  over  that  expected  for 
a  Mach  3.  1  normal  shock  wave  can  be  explained  by  the  momentum  loss  on 
the  fuel  injection  wedge. 


3.1.2  Pressure  Effects  Aft  of  the  Normal  Shock  Wave 


Several  pressure  profiles  were  taken  in  the  region  behind  the  normal 
shock  wave  with  the  preheater  on  but  without  test  fuel  injection.  Figure  8 
shows  the  data  from  an  axial  static  pressure  and  a  pitot  total  pressure 
traverse.  The  static  pressure  at  the  shock  wave  was  calculated  from  the 
theoretical  normal  shock  static  pressure  rise  at  the  centerline;  based  on 
the  measured  indicated  total  pressure  and  the  isentropic  static  pressure  at 
a  free- stream  Mach  number  of  3.1,  the  static  pressure  falls  rapidly  from 
the  shock  wave  aft  and  the  indicated  total  pressure  falls  more  slowly. 

As  may  be  seen  from  the  data  (Fig.  8),  the  static  pressure  is  almost 
independent  of  total  temperature.  This  is  not  true  for  the  total  pressure. 
As  the  total  temperature  is  increased,  the  indicated  total  pressure  ratio 
increases  (Fig.  9).  This  increase  results  from  a  greater  spreading  of  the 
wake  from  the  fuel  injection  wedge  at  higher  temperatures  [see  fuel  dis¬ 
tributions  (Fig.  14)]  which  spreads  the  wedge  loss  over  a  greater  area  and 
reduces  the  maximum  loss. 

Pressure  profiles  taken  in  the  vertical  and  horizontal  direction  re¬ 
sulted  in  curves  of  the  type  seen  in  Fig.  10.  The  horizontal  profiles  are 
nearly  flat  in  the  center.  The  vertical  profiles  are  strongly  curved  indi¬ 
cating  both  the  disturbance  from  the  wedge  and  the  transition  in  Mach 
number  as  the  probe  crosses  the  slip  lines  from  the  shock  intersection. 
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From  the  static  and  indicated  total  pressure  data,  the  centerline 
Mach  number  and  total  pressure  may  be  calculated.  These  parameters 
are  plotted  in  Figs.  8  and  9.  The  Mach  number  increases  from  the 
normal  shock  value  just  behind  the  shock  wave  to  a  supersonic  value  at 
less  than  three  inches  downstream  of  the  shock  wave  (Fig.  8).  The 
Mach  number  at  three  inches  aft  also  increases  with  an  increase  in  total 
temperature  (Fig.  9). 

In  SCT  when  a  combustible  mixture  is  heated  by  passing  through 
the  normal  shock  wave,  it  immediately  enters  a  flow  field  where  the 
velocity  is  increasing  and  the  static  temperature  is  decreasing.  There¬ 
fore,  any  pressure  losses  from  combustion  will  depend  greatly  on  the 
axial  position  at  which  the  combustion  occurs.  Also,  the  static  tem¬ 
perature  drop  from  the  flow  acceleration  will  affect  the  temperature 
rise  from  combustion  and  may  greatly  change  the  course  of  the  reaction. 


3.1.3  Total  Temperature  Profiles 

Total  temperature  profiles  were  made  at  various  locations  in  the 
test  section.  The  temperatures  were  determined  from  the  analysis  of 
the  combustion  gases  from  the  hydrogen-fired  preheater  using  100- 
percent  combustion  efficiency  (see  Procedure).  No  free  hydrogen  was 
ever  detected  in  the  test  section  unless  test  fuel  was  being  introduced. 
Typical  temperature  profiles  may  be  seen  in  Fig.  11.  The  profiles  are 
relatively  flat  over  the  area  of  interest  at  any  given  axial  station.  There 
is,  however,  a  gradual  axial  drop  in  temperature  toward  the  rear  which 
results  from  mixing  of  the  heated  core  with  surrounding  cooler  air. 

As  the  total  temperature  is  increased,  the  oxygen  content  of  the  air 
decreases  and  the  water  vapor  increases.  Figure  12  shows  the  ideal 
relationship  between  the  total  temperature  and  the  gas  composition.  This 
produces  an  added  complication  since  the  physical  and  chemical  properties 
of  the  preheater  gases  are  a  very  complicated  function  of  temperature. 
Figure  13  shows  how  some  of  the  more  useful  of  these  properties  vary. 

For  the  preheater  gases  without  test  fuel,  Tt^  is  the  total  temperature 
at  the  sampling  point.  When  test  fuel  is  added,  Tt^  is  the  total  temper¬ 
ature  of  the  mixture  with  «  =  0.  it  is  the  final  temperature  of  the  combus¬ 
tion  products  at  the  same  location. 

The  addition  of  cold  test  fuel  changes  the  temperature  and  composition 
still  further.  This  effect  will  be  discussed  in  the  section  on  combustion 
efficiency. 
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3.2  FUEL  DISTRIBUTION 

Fuel  injected  from  a  small  hole  in  the  downstream  side  of  the  fuel 
injector  wedge  (Fig.  3)  proceeds  downstream.  While  traveling,  it 
mixes  with  the  surrounding  air.  The  rate  of  mixing  is  influenced  by 
difference  in  velocity  between  fuel  and  air  and  by  temperature  levels. 
After  traveling  about  six  inches  in  a  Mach  3  flow,  the  mixture  passes 
through  a  shock  where  reaction  may  or  may  not  occur  depending  on  the 
gas  temperature. 

Fuel  concentration  was  calculated  from  analysis  of  oxygen  and  hydro¬ 
gen  content  by  the  method  described  in  Appendix  A.  Data  with  test  fuel 
combustion  will  give  the  hydrogen  equivalence  ratio  (y^  +72  +  3^).  Since 
the  same  preheater  level  without  test  fuel  will  give  y^,  the  test  fuel  will 
be  the  difference  of  the  two  calculated  equivalence  ratios.  The  resulting 
fuel  concentration  profiles  were  plotted  in  Fig.  14  and  closely  resemble 
a  normal  distribution  curve  of  the  form  K(e)  -i>z\  where  K  and  b  are 
coefficients  and  Z  is  a  tunnel  coordinate  distance  perpendicular  to  the 
flow  in  the  horizontal  plane.  Figure  14  shows  a  marked  increase  in  rate 
of  mixing  of  the  fuel  as  preheater  temperature  is  increased.  A  compar¬ 
ison  of  fuel  profile  peaks  at  the  fore  and  aft  position  in  Fig.  14  also  shows 
a  decrease  in  fuel  concentration  with  distance  from  the  normal  shock. 


3.3  COMBUSTION  EFFECTS 

One  of  the  primary  purposes  of  these  experiments  was  to  show  that 
heat  release  actually  occurred  in  the  area  aft  of  the  shock,  and,  if  pos¬ 
sible,  to  what  extent.  This  was  investigated  by  three  methods;  temper¬ 
ature  measurements,  pressure  measurements,  and  radiation  emission. 
They  will  be  discussed  separately. 

3.3.1  Effect  of  Temperature  on  Combustion 


One  of  the  overall  objectives  of  this  study  was  to  determine  the  rate 
of  reaction  of  hydrogen-air  mixtures.  As  an  initial  part  of  this  investi¬ 
gation,  the  combustion  efficiency  was  measured  at  various  locations  in 
the  combustion  zone.  These  data  are  correlated  with  position  relative  to 
the  shock  wave  and  with  the  inlet  temperature.  An  attempt  was  also  made 
to  determine  a  reaction  rate  as  a  function  of  the  combustion  mixture 
temperature. 

3.  3.  1.  1  Combustion  Efficiency.  The  combustion  efficiency,  f ,  is 
defined  in  this  report  as  the  ratio  of  the  test  fuel  burned  to  the  total  test 
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fuel  present  that  could  burn.  As  was  previously  mentioned  the  preheater 
operates  at  100-percent  combustion  efficiency,  and  the  preheater  hydrogen 
is  not  counted  in  calculating  the  test  fuel  combustion  efficiency. 

Accuracy  with  which  the  efficiency  and  temperature  can  be  meas¬ 
ured  depends  on  two  factors:  (1)  the  accuracy  with  which  the  composition 
of  the  sampled  gas  can  be  determined  (see  Appendix  A)  and  (2)  the  rate 
at  which  the  reaction  is  quenched  as  it  enters  the  probe.  The  rate  of 
quench  is  discussed  in  greater  length  in  Appendix  B.  However,  it  should 
be  noted  that  free  radical  recombinations  cannot  be  quenched;  that  is, 
any  OH,  H,  or  O  will  recombine  to  form  H2O,  H2,  and  02.  Any  energy 
released  during  these  recombinations  will  increase  the  temperature  of 
the  gas  after  it  has  entered  the  probe.  The  recombinations  can  be 
ignored  if  the  efficiency  is  considered  as  the  hydrogen  consumed.  Calcu¬ 
lated  temperatures  and  energy  release  will  be  higher  than  the  free- stream 
values  by  the  dissociation  energy  of  the  free  radicals  present,  even  if  it 
is  assumed  that  the  dissociation  and  branching  reactions  are  completely 
quenched. 

The  relationship  between  e  three  inches  downstream  of  the  normal 
shock  and  total  temperature  (Tiq)  is  shown  in  Fig.  15.  Figures  15a,  b, 
and  c  show  how  e  varies  along  the  horizontal  centerline  for  several  inlet 
temperatures.  The  gas  temperature  from  the  preheater  is  nearly  uni¬ 
form  along  the  horizontal  centerline  when  no  test  fuel  is  added;  but  in 
calculating  a  correction  for  the  cold  test  fuel  it  was  found  that  there  is  a 
reduction  in  total  temperature  which  is  proportional  to  the  total  test  fuel 
added.  Figure  16  shows  «  as  a  function  of  total  temperature  (Tt^)  cor¬ 
rected  for  the  cooling  effect  of  the  hydrogen  at  an  axial  position  three 
inches  aft  of  the  normal  shock  wave.  At  this  station  for  a  plenum  total 
pressure  of  50  psia  the  efficiency  is,  within  the  experimental  error,  a 
function  of  inlet  total  temperature  only,  with  no  noticeable  effect  of  fuel 
concentration.  The  static  temperature  of  the  gas  just  behind  the  shock 
wave  may  be  calculated  by  multiplying  Tt^ ,  corrected  for  the  total  tem¬ 
perature  rise  along  the  axis  (Fig.  11c),  by  the  static  to  total  temperature 
ratio  at  the  shock  wave  (From  M  in  Fig.  8).  A  resulting  inlet  static 
temperature  of  about  2300°R  is  required  to  produce  50-percent  com¬ 
bustion  efficiency  in  three  inches  and  about  1950°R  to  produce  five-percent 
combustion  in  the  same  distance. 

In  a  system  such  as  this  where  the  flow  accelerates  rapidly  behind 
the  shock  wave,  the  combustion  efficiency  rises  rapidly  at  first,  and  then 
as  the  static  temperature  reaches  a  maximum,  its  rate  of  rise  decreases 
(Fig.  17).  The  static  temperature  rises  rapidly  at  first  because  the  total 
temperature  is  increasing.  However,  at  the  same  time  the  static  to  total 
temperature  ratio  is  decreasing  because  of  the  flow  acceleration.  As  the 
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reaction  nears  completion  this  effect  becomes  predominant,  and  the 
static  temperature  begins  to  decrease.  Since  the  static  temperature 
was  calculated  from  Mach  number  and  total  temperature,  the  accu¬ 
racy  of  this  parameter  depends  on  the  accuracy  of  both  the  total  tem¬ 
perature  and  the  indicated  static  and  total  pressure  measurements. 

3.  3.  1.  2  Reaction  Kinetics.  The  relation  between  combustion  effi¬ 
ciency  and  distance  can  be  transformed  to  a  relationship  between  combus¬ 
tion  efficiency  and  time  by  calculating  the  velocity- distance  relationship 
from  experimental  values  of  Pg.  Pt^  and  Tt  and  graphically  integrating 

1/V  vs  X  where  /  i/v  dX  =  r{X).  Figure  18  shows  the  f  plotted  as  a 
0 

function  of  r.  The  slope  of  this  curve  is  dt/dr  which  equals  -  - - 

(Hj)j  dr 

A  global  equation  for  the  rate  of  consumption  of  hydrogen  can  be  written  as 


d  [Hj 
d 


[H,] 


[OJ. 


where  A  is  a  constant,  E  is  the  global  activation  energy,  and  m  and  n 
are  constants  which  determine  the  overall  order  of  the  reaction. 

If  the  proper  values  for  m  and  n  are  known  and  one  relationship  will 
describe  the  whole  reaction,  a  plot  of  log  Ae  E/RoT 

1/T  will  be  a 

straight  line,  where  the  slope  is  E/Rq  and  the  intercept  at  E/RoT  =  1 
is  A  (Fig.  19).  The  data  shown  in  this  figure  came  from  two  runs  at 
different  inlet  temperatures.  The  activation  energies,  as  represented 
by  the  slope  of  the  curves,  were  calculated  from  the  two  sets  of  data  and 
are  nearly  the  same  above  a  static  temperature  of  2800°R.  The  change 
in  slope  below  2800°R  could  result  because  the  global  equation,  as  written, 
does  not  describe  the  reaction  over  a  broad  temperature  range  or  because 
of  errors  in  the  data  at  the  lov/er  temperature.  The  difference  in  the  inter 
cepts  of  the  two  curves  also  shows  that  the  equation  does  not  describe  the 
reaction  completely.  The  calculated  activation  energy  is  close  to  that  for 
the  reaction  H2  +  O  -*•  OH  +  H  which  has  been  proposed  as  the  rate  con¬ 
trolling  step  in  the  hydrogen- oxygen  combustion  (Ref.  4).  However,  there 
are  an  insufficient  number  of  data  points  to  establish  a  definite  value  for 
the  constants  in  a  global  equation. 

Between  0.  5  in.  and  3  in.  aft  of  the  shock  wave,  the  change  in  com¬ 
bustion  efficiency  is  very  small  (Fig.  17).  Reaction  rates  in  this  range 
are  much  lower  and  do  not  fit  the  same  global  equation  as  they  do  in  the 
early  part  of  the  reaction.  It  is  possible  that  the  rate  of  change  of  hydro¬ 
gen  concentration  is  affected  by  mixing  or  that  an  apparent  change  results 
from  a  failure  of  the  probe  to  follow  a  streamline. 
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3.3.2  Effect  of  Temperature  Rise  on  Total  Pressure 

As  a  further  attempt  to  verify  heat  release,  total  pressure  was 
measured  at  various  inlet  total  temperature  levels  and  locations  in  the 
zone  of  interest  downstream  of  the  normal  shock,  both  with  and  without 
test  fuel.  Also,  a  simplified  theoretical  correlation  of  this  data  was 
attempted.  For  the  case  of  no  friction,  mass  addition,  or  molecular 
weight  change  (Ref.  6),  the  following  relation  may  be  written  between 
total  pressure  and  total  temperature 

dpt  yM^  dT^ 


which,  when  integrated  for  the  condition  of  constant  Mach  number  and  y. 
becomes 


Although  total  pressure  ratio  is  a  function  not  only  of  total  temper¬ 
ature  ratio,  but  also  of  Mach  number  and  y»  Eq.  (2)  may  be  used  for  a 
practical  case  where  total  temperature  increases  are  small,  and  y, 

Mach  number,  mass,  and  molecular  weight  changes  may  be  considered 
to  be  insignificant.  With  this  equation  in  mind,  the  data  on  pressure 
measurement  will  be  examined. 

First,  Fig.  15  shows  there  are  two  pressure  profile  curves  for  each 
inlet  gas  temperature  level  (Tt^),  one  for  the  condition  of  preheater  only 
and  the  other  for  preheater  plus  test  fuel.  Note  also  that  the  difference 
between  the  two  curves  increases  with  increased  inlet  temperature  level. 
Here  a  total  pressure  loss  is  apparent  when  combustion  occurs;  however, 
total  pressure  loss  may  result  from  reduced  molecular  weight  also. 

Then  in  Fig.  20,  the  total  pressure  drop  caused  by  unburned  fuel  is 
shown,  and  in  Fig.  21  the  effect  of  combustion  is  demonstrated  by  adding 
in  a  large  group  of  data  at  various  degrees  of  reaction  completion.  Now 
a  family  of  lines  can  be  drawn,  showing  the  relation  between  heat  release, 
fuel  concentration,  and  total  pressure  loss.  If  the  effect  of  pressure  drop 
caused  by  unburned  fuel  is  removed  by  using  the  data  of  Fig.  20,  the  curve 
of  Fig.  22  may  be  drawn. 

Thus  far,  only  experimental  data  have  been  discussed,  and  it  can  be 
seen  that  the  scattering  of  these  results  are  well  within  the  estimated 
experimental  accuracy,  as  defined  in  Appendix  A.  In  calculating  the  points 
for  these  curves,  it  was  necessary  to  use  total  and  static  pressure  measure¬ 
ments  and  total  temperature,  y,  and  molecular  weight  from  the  gas  analysis. 
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Now,  the  question  arises,  can  the  raeasurement  of  total  pressure 
only  be  used  to  predict  the  quantity  of  heat  release.  The  answer  is  no 
because  enough  additional  information  must  be  known  to  describe  the 
properties  of  the  flowing  gas,  such  as  Mach  number,  y,  and  molecular 
weight.  However,  enough  experimental  information  was  obtained  so  that 
the  pressure  drop  caused  by  heat  addition  could  be  separated  from  that 
caused  by  fuel,  and  with  what  is  considered  good  agreement. 

Equation  (2)  was  used  with  a  stepwise  calculation.  Temperature 
increments  of  100°R  were  used.  The  initial  Mach  number  conditions 
were  calculated  just  aft  of  the  normal  shock,  based  on  the  plenum  and 
test  section  pressures,  analysis  of  the  wedge  effect  on  the  flow  (Sec¬ 
tion  3.  1.  1),  and  the  fuel  concentration.  The  axial  static  pressure  pro¬ 
file  was  experimentally  found  to  be  independent  of  heat  release  at 
distances  greater  than  one  inch  aft  of  the  normal  shock  wave.  Just  behind 
the  shock  wave  the  static  pressure  was  calculated  from  the  free- stream 
static  pressure  at  Mach  3.  1  and  the  measured  total  pressure.  Using  this 
value  and  the  experimental  data,  the  pressure  distribution  curve  of 
Fig.  23  is  shown.  The  specific  heat  ratio  (y)  was  calculated  from  the  gas 
analysis  and  temperature  determined  from  the  gas  analysis.  Thus,  an 
incremental  total  pressure  drop  was  calculated,  resulting  in  a  new  total 
pressure  for  the  next  increment.  In  this  manner,  the  calculation  was 
continued  until  the  temperature  rise  tapered  off  to  zero  at  a  distance  about 
0.  6  in.  downstream  of  the  shock.  A  comparison  between  the  calculated 
total  pressure  loss  and  the  experimental  data,  for  the  condition  of  test 
fuel  injected  and  preheater  on,  is  given  in  Fig.  24  at  two  different  inlet 
temperatures.  This  correlation  is  considered  good  because  a  Mach  num¬ 
ber  error  of  0.  05  used  in  the  stepwise  calculation  would  result  in  an  error 
of  about  0.  2  percent  per  step  and  an  overall  error  of  about  3  percent. 

The  uncertainties  in  the  static  pressure  measurement  could  easily  cause 
this  difference. 

These  experimental  pressure  data  were  corrected  for  unburned  fuel. 

It  was  assumed  that  the  difference  between  the  fuel-on  and  fuel- off  total 
pressure  just  aft  of  the  normal  shock  (X  =  0)  was  caused  by  the  presence 
of  unburned  fuel,  and  that  as  the  fuel  burned,  this  effect  was  reduced  pro¬ 
portionally  to  the  concentration  of  fuel.  Although  water  vapor  was  pres¬ 
ent,  its  overall  effect  was  small  compared  to  that  of  hydrogen,  and  there¬ 
fore  it  was  neglected.  The  resulting  equations  are: 


_  Pressure  loss  ratio  caused  by  fuel  (3) 

just  aft  of  the  normal  shock 

X  =  0 


(y.  -  (^).. 
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and 


[(y,  +  a^)  -  y,^.]  X 

(y^  + 


so  that 


ID  e  a  s 


eat  add, 
only 


1  - 


X 


(4) 


(5) 


The  terms  in  the  preceding  equations  are  located  on  the  following 
schematic  diagram: 


The  correlation  was  calculated  only  to  about  0.  6  in.  aft  of  the  normal 
shock  since  (1)  most  of  the  temperature  rise  occurs  within  this  distance  and 
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(2)  it  was  surmised  that  mixing  losses  were  becoming  pronounced  at  this 
point,  which  would  have  made  further  calculations  irrelevant  since  the 
equations  used  did  not  account  for  mixing. 

From  the  foregoing  data  and  analysis,  these  conclusions  may  be 
stated; 

1.  There  is  a  well  defined  pressure  loss  effect  that  can  be  attrib¬ 
uted  to  (1)  fuel  injection  losses  when  hydrogen  is  added  and  to 
(2)  heat  addition  from  combustion. 

2.  This  effect  may  be  correlated  with  simple  theory,  even  though 
the  flow  field  is  changing,  if  sufficient  information  on  the  flow 
field  is  at  hand. 

3.3.3  Test  Section  Emission 


During  the  investigation  of  hydrogen  combustion  initiated  by  a  shock 
wave,  pronounced  emission  was  visible  in  the  region  behind  the  shock 
wave . 

3.  3.  3.  1  Normal  Shock  Wave.  An  emission  photograph  downstream 
of  a  normal  shock  wave  was  superimposed  on  a  schlieren  photograph  and 
is  shown  in  Fig.  25.  This  type  of  emission  has  been  used  to  calculate  an 
ignition  delay  time  (Refs.  1,  2,  and  4).  However,  in  the  SCT,  although 
emission  from  the  test  section  begins  at  an  inlet  total  temperature  of 
about  1700®R,  the  combustion  efficiency  does  not  exceed  five  percent  unless 
the  total  temperature  is  above  2000®R. 

In  an  attempt  to  determine  the  source  of  this  radiation,  a  spectrogram 
(covering  a  wavelength  from  about  3000  to  6000  A)  was  taken  of  the 
brightest  part  of  the  emission  zone.  A  quartz  prism  spectrophotometer 
with  a  photomultiplier  detector  was  used. 

O 

The  only  radiation  which  could  be  detected  was  emitted  around  3050  A 
which  corresponds  to  the  band  head  of  the  0,  0  OH  band  at  3064  A,  and  con¬ 
tinuous  radiation  which,  when  analyzed,  had  the  spectral  distribution  of  a 
black  body  at  3400°R.  The  most  logical  explanation  for  the  continuous 
radiation  is  that  it  comes  from  particulate  matter  in  the  stream.  This  is 
not  to  say  that  the  particles  are  at  3400‘’R,  since  the  emissivity  of  small 
particles  will  vary  with  wavelength  and  the  relationship  between  true  tem¬ 
perature  and  apparent  black  body  temperature  is  quite  complicated.  The 
particles  apparently  were  heated  by  the  hydrogen  combustion,  since  at  inlet 
total  temperatures  up  to  2800°R  there  is  no  glow  without  hydrogen.  Above 
2800°R  there  is  an  overall  glow  in  the  test  section  which  becomes  brighter 
as  the  temperature  is  increased  and  which  is  brighter  behind  the  shock  waves 
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where  the  static  temperature  is  higher  (Fig.  26).  The  fact  that  the  loca¬ 
tion  where  the  emission  brightens  is  coincidental  with  the  location  of  the 
shock  waves  would  imply  that  at  least  some  of  the  particles  are  very  small 
(<<  1  micron)  since  they  must  be  heated  appreciably  in  less  than  one  micro¬ 
second.  Any  longer  delay  than  this  could  be  detected  as  a  gap  between  the 
shock  waves  as  shown  by  the  schlieren  system  and  the  shock  waves  as  out¬ 
lined  by  the  emission  increase.  Other  particles  must  be  substantially 
larger  since  there  is  emission  ahead  of  the  shock  wave  where  the  static 
temperature  is  about  1000°R,  and  the  particles  radiating  here  must  have 
retained  their  heat  from  the  preheater.  Within  the  limits  of  knowledge  of 
the  temperatures  of  the  particles,  a  particle  size  distribution  from  10  to 
0.  1  micron  would  probably  satisfy  these  conditions.  Gas  stream  analysis 
has  shown  small  quantities  of  iron  and  aluminum  oxides  to  be  present  in 
the  supply  air.  This  contamination  probably  originates  from  scale  in  the 
pipe  and  from  the  aluminum  oxide  air  drier  and  could  account  for  the  ob¬ 
served  effect. 

It  is  possible  that  the  emission  which  is  seen  could  result  from  sur¬ 
face  reactions  which  heat  the  particles  without  appreciably  changing  the 
gas  temperature  or  the  composition  of  the  bulk  of  the  stream.  It  is  not 
known  whether  particles  heated  by  surface  reactions  would  affect  the 
overall  rate  of  the  reaction,  but  it  is  possible  that  these  reactions  could 
cause  an  increase  in  the  concentration  of  free  radicals  which  could  change 
the  rate  of  the  branching  reactions.  Also  later  in  the  reaction,  the 
particles  could  provide  a  surface  for  recombination  reactions.  Both  of 
these  effects  would  cause  an  increase  in  the  overall  rate  of  conversion 
of  hydrogen  and  oxygen  to  water  vapor. 

3.  3.  3.  2  Oblique  Shock  Wave.  The  emission  seen  behind  the  shock 
wave  when  a  normal  shock  configuration  was  used  was  also  seen  with  the 
oblique  shock  wave  configuration.  Figure  27  shows  emission  photographs 
of  the  test  section  at  an  inlet  total  temperature  of  over  3000°R  with  nitro¬ 
gen  and  then  with  hydrogen  flowing  through  the  fuel  injector.  Both  photo¬ 
graphs  were  processed  in  the  same  manner  and  show  the  relative  intensity 
and  location  of  the  emission  from  the  hot  preheater  gases  and  from  the 
hydrogen  reaction.  The  calculated  static  temperature  immediately  behind 
the  oblique  wave  is  the  same  as  that  behind  a  normal  shock  wave  with  a 
total  inlet  temperature  of  about  2300°R.  Under  these  conditions  a  reaction 
should  proceed  in  the  main  body  of  the  gas  on  the  order  of  40  percent  of 
the  hydrogen  burned  in  the  first  inch  behind  the  oblique  shock  Wave.  This 
conclusion  is  based  on  the  effect  of  inlet  temperature  on  combustion 
efficiency  (Fig.  16)  and  on  the  spacial  distribution  of  the  combustion  zone 
(Fig.  17). 

Several  attempts  were  made  to  determine  if  there  were  a  Mach  num¬ 
ber  change  resulting  from  heat  addition  behind  an  oblique  wave  by  looking 
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for  wave  angle  increases  on  small  wedges  mounted  in  the  emission  zone. 
Results  of  the  data  were  inconclusive  because  of  (1)  aerodynamic  effects 
which  prevented  the  shock  waves  from  being  straight  and  clear-cut  and 
(2)  reduced  density  at  high  temperature  which  reduced  the  visible  shock 
wave  intensity  as  seen  in  the  schlieren  system.  No  data  on  composition 
changes  were  obtained  because  the  gas  sampling  system  was  not  in  use 
at  the  time. 


4.0  CONCLUDING  REMARKS 


A  series  of  tests  were  conducted  in  the  Supersonic  Combustion  Tun¬ 
nel  to  determine  the  efficiency  of  the  combustion  in  a  shock- induced 
combustion  wave,  the  effect  of  the  combustible  mixture  temperature  on 
the  combustion  efficiency,  and  the  relationship  between  the  total  pres¬ 
sure  loss  and  the  amount  of  heat  release. 

An  aerodynamic  analysis  of  the  Mach  reflected  normal  shock  system 
showed  a  static  pressure  gradient  behind  the  normal  shock  wave  which 
caused  the  flow  to  be  re-accelerated  to  about  Mach  1.  6  three  inches 
downstream  of  the  normal  wave.  This  pressure  gradient  was  independ¬ 
ent  of  total  temperature  level  and  was  unaffected  by  the  presence  of 
combustion  behind  the  shock  wave. 

Hydrogen  for  shock- induced  combustion  was  introduced  from  the 
center  of  a  wedge-shaped  strut  located  in  the  supersonic  flow  about  six 
inches  upstream  of  the  normal  shock  wave.  The  fuel  spread  so  that  the 
entire  height  of  the  normal  shock  wave  contained  a  combustible  mixture. 
This  spreading  increased  further  downstream  and  also  increased  as  the 
total  temperature  was  raised. 

When  test  fuel  was  introduced,  the  combustion  efficiency  rose  rapidly 
at  first  and  then  at  a  lower  rate.  The  decrease  in  rate  resulted  from  a 
reduction  in  the  quantity  of  hydrogen  available  to,  burn  and  from  the  re¬ 
duction  in  static  temperature  caused  by  the  expansion  of  the  gas.  Three 
inches  downstream  of  the  shock  wave  the  combustion  efficiency  was  a 
function  of  inlet  total  temperature  and  independent  of  the  initial  hydrogen 
concentration  within  the  experimental  accuracy  of  the  data.  Initial  com¬ 
bustion  occurred  at  a  total  temperature  of  about  1800°R,  increased  to 
50  percent  at  about  2300°R  and  to  90  percent  at  about  SOOO^R. 

The  effect  of  temperature  rise  on  total  pressure  loss  was  measured 
for  several  conditions  of  inlet  temperature  (Tt„)  and  position  downstream 
of  the  normal  shock  and  was  correlated  with  fuel  concentration  and  com¬ 
bustion  efficiency.  In  addition,  a  correlation  of  pressure  loss  and 
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temperature  rise  was  obtained  by  calculating  the  pressure  loss  for  an 
incremental  temperature  rise  with  a  simplified  equation  and,  using  the 
known  factors  about  the  flow  field,  then  comparing  the  resulting  values 
with  the  experimental  data.  Correlation  between  the  calculated  and  ex¬ 
perimental  values  was  shown  in  the  region  just  aft  of  the  normal  shock 
where  the  temperature  rise  occurred. 

Emission  from  the  combustion  zone  was  detected  at  inlet  total  tem¬ 
peratures  as  low  as  1700°R  and  increased  in  brightness  as  the  temper¬ 
ature  was  raised.  Spectral  analysis  showed  this  emission  to  consist  of 
some  radiation  about  the  wavelength  of  the  OH  band  emission  and  con¬ 
tinuous  radiation  which  increased  from  short  to  long  wavelengths.  It 
is  possible  that  the  continuous  emission  came  from  particles  of  iron  and 
aluminum  oxides  heated  to  incandescence  by  surface  chemical  reactions. 

Tests  were  made  with  a  single  wedge  in  an  attempt  to  get  oblique 
shock-induced  combustion.  Emission  was  seen  behind  the  shock  wave, 
but  since  no  chemical  analysis  of  the  gas  was  taken  behind  the  shock 
wave  no  quantitative  evidence  of  combustion  was  obtained. 
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APPENDIX  A 

CALCULATION  PROCEDURES 


CHEMICAL  REACTIONS 

The  chemical  reactions  of  interest  are  located  physically  in  two 
places:  in  the  preheater  where  1460°R  air  is  heated  to  a  higher  temper¬ 
ature  by  combustion  of  hydrogen  and  in  the  test  section  where  fuel  is 
injected  and  reaction  phenomena  are  observed. 

Since  all  data  to  date  indicate  no  residual  hydrogen  from  the  pre¬ 
heater  reaction,  combustion  efficiency  there  was  considered  to  be 
100  percent.  Combustion  efficiency  in  the  test  zone  is  based  on  the 
ratio  of  fuel  quantity  that  has  disappeared  in  a  reaction  process  as  com¬ 
pared  to  the  total  fuel  that  could  theoretically  react.  Hence,  if  the  fuel 
concentration  is  greater  than  stoichiometric  (equivalence  ratio  >1),  the 
quantity  of  fuel  in  excess  of  stoichiometric  is  automatically  excluded 
from  the  efficiency  calculation.  In  order  to  accomplish  this  calculation, 
it  was  found  necessary  to  calculate  overall  combustion  efficiency  of  pre¬ 
heater  and  test  fuel,  insert  the  known  value  of  preheater  fuel  concentra¬ 
tion  at  100  percent  combustion  efficiency,  then  calculate  the  resulting 
combustion  efficiency  of  the  test  fuel. 

The  equilibrium  reaction  equation  of  H2"air  after  quenching  and 
cooling,  for  various  efficiencies,  is  as  follows; 

2(y,  +  y,  +  a,)  H,  +  0,  +  3.77  N,  - 

[1  -  (y^  +  y^)r}]  Oj  +  2(y^  +  y^)  T]  H^O 

+  2(y.  +  y^)  (1  -  r?)  H,  +  2a,H2  +  3.77  N, 

where  1  (overall  combustion  efficiency)  is  defined  as: 

amount  of  fuel  burned 

71  =  —  .  . .  . .  -  - - 

amount  of  fuel  that  could  be  burned  with  available  O2 
«  =  combustion  efficiency  for  test  zone  fuel 

yi  =  equivalence  ratio  of  fuel  from  preheater 

yj  =  equivalence  ratio  of  test  fuel 

aj  =  equivalence  ratio  of  test  fuel  greater  than  E.R.  =  1.0 
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The  completely  burned  gas  from  the  preheater  may  be  analyzed  and 
equivalence  ratio  (E.K.)  calculated  from  Eq.  (2)  where  H2%  =  0  and 

■q  =  100%: 


1  -  4.77  (Oj) 

~  1  -  Oj 

where 

(Oj)  =  percent  of  oxygen  in  the  gas  by  volume, 
on  a  dry  basis 


{A-2) 


Once  Vi  is  established  at  a  fixed  preheater  level,  test  fuel  can  be 
injected.  Since  test  fuel  can  exist  from  0  to  100  percent  burned,  both 
oxygen  and  hydrogen  analyses  are  required  to  determine  degree  of  the 
reaction  completion.  Here,  it  is  assumed  that  the  reaction  is  quenched 
a  short  interval  of  time  after  entering  the  sampling  probe  (see  Appen¬ 
dix  B),  and  that  free  radicals  have  recombined  with  their  proper 
associates  during  the  cooling  process.  With  these  assumptions,  the 
following  equation  may  be  obtained  from  Eq.  (A-1): 


(yi  +  72  +  aj)  = 


1  -  4.77  (O2)  +  0.88  (H2) 
1  -  [  (H,)  +  (O2)  1 


(A-3) 


where  (H2)  -  percent  of  hydrogen  by  volume,  dry. 

Since  is  not  known,  (y^  +  32)  may  be  calculated.  Since  aj  can  be 
greater  than  zero  only  if  (y^  +  Xj  +  82)  >  1 ,  then  a2  and  y^  may  be  evaluated. 


From  Eq.  (A-1)  may  also  be  derived  the  equation  for  ??: 

_ 1  -  4.77  (O2)  -  (H2) _ 

^  1  -  4.77  (O2)  +  0.88  (Hj)  +  [  (H2)  +  (O2)  -  1  ]  (aj) 


(A-4) 


Combustion  efficiency  for  the  test  fuel  only  may  then  be  defined  as; 

V  ivi  +  72)  “  Xi  (A- 5) 

yj 


GAS  TEMPERATURE  CALCULATION 

In  calculating  gas  temperature,  dissociation  was  assumed  to  be 
negligible.  Since  dissociation  is  not  appreciable  until  temperatures 
above  4000‘’R  and  near  or  above  stoichiometric  ratios  are  reached,  it 
was  neglected  for  the  data  of  this  report. 
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The  final  gas  temperature  was  obtained  from  the  following  equations: 

Enthalpy  Balance 


(™fH,  hair  +  ^  Qr  )  = 

(A-6) 

hH,  +  mfo^  ho,  +  hN,  + 

On  the  basis  of  one  mole  of  0,  the  mole  function  of  the  reactants  and 
products  are; 

Reactants  Products 

"•fH,  =  2  +  72  +  =  2  (y^  +  )  (1  -  7? )  +  2  a, 

""fair  =  4.77  =  1  _  (y^  +  yjjy 

“fN,  =  3*77 

"’fH,0  =  2^  (7i  +  7^) 


Enthalpy  equations  for  the  products  in  the  form  h  =  a  +  bTt  +  cTt^ 
were  developed  for  H, ,  0,,  iN,,  and  H,0  from  1200  to  4000°R  using  enthalpy 
data  from  Ref.  13: 

hn,  =  -3184  +  6.2596  Tt  +  0.2683  x  10"'  Tt'  Btu/lb  mole  ( +  10) *  (A-7) 

ho,  =  -2818  +  6.1636  Tt  +  0.4786  x  10“'  Tt'  Btu/lb  mole  (+  300)  (A-8) 

hN^  =  -3941  +  6.9532  Tt  +  0.2411  x  10“'  Tt'  Btu/lb  mole  (+40)  (A- 9) 

hHjO  =  -4277  +  7.441  Tt  +  0.682  x  10“'  Tt'  Btu/lb  mole  (±60)  (A-10) 

In  these  equations,  h  =  0  at  524°R. 

Similar  equations  for  the  reactants  were  developed: 

hR,  =  -3686  +  6.984  Ttjj^  (±10)*  Tt  =  600  -  1000°R  (A- 11) 

hair  =  -4348  +  7.617  Tt^.,  (±20)  Tt  =  1300  -  1600°R  (A- 12) 

where  h  =  0  at  524°R  and  Qj^  =  104100  Btu/lb  mole  at  524°R.  (Ref.  14) 


*The  equation  fits  the  tables  (Ref.  13)  to  this  tolerance  expressed 
as  Btu/lb  mole. 
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These  equations  may  be  combined  to  form  the  quadratic  equation; 
[  1.3876  +  0.5366  (y^  0.3488  r/  (y^  +  y^)  ]  10  ^  Tt^ 

+  [  32.3772  +  12.5192  (y^  +  +  a,)  -  3.8008  r?  (y^  +  y^)  ]  Tt  (A- 13 

-  [-3063  -  998  (y^  +  Xj  +  a^)  +  207600  rj  {y^  +  y^) 

+  36.333  Tt^.^  +  13.968  (y^  +  y^  +  a,)  ]  =  0 


These  equations  for  a  Hj  -  air  reaction  mixture  are  accurate  to  ±5‘’R. 

Three  gas  temperatures  may  be  calculated: 

1,  Preheater  only  (Tt^)  where  and  aj  are  0  and  rj  =  1.0. 

2,  Preheater  plus  test  fuel  without  test  fuel  combustion  (.Tt^)  where 
V  (y,  +  yj  =  y^. 

3,  Preheater  plus  test  fuel  with  test  fuel  combustion  (Tt)  as  written 
in  Eq.  (A- 13). 

Thus,  the  properties  of  the  reacting  gas  that  can  be  determined  from 
a  Hj  -  O2  analysis  on  a  dry  basis,  assuming  a  quenched  reaction,  are  as 
follows: 

1.  Fuel-air  ratio  (partially  burned  or  unburned). 

2.  Combustion  efficiency. 

3.  Gas  temperature  from  fuel-air  ratio  and  combustion  efficiency. 

ERROR  ANALYSIS 

Error  analysis  was  based  on  an  estimated  95 -percent  probability 
that  the  measurements  lay  inside  a  certain  interval.  For  measured  data, 
the  interval  was  taken  as  that  within  which  repeatable  observations  could 
be  made.  For  calculations  the  plus  or  minus  increment  (A)  was  calcu¬ 
lated  from  the  data  a’s  .  A  comparison  of  Table  A-1  with  Figs.  16,  19, 

20,  and  21  will  verify  the  validity  of  the  estimated  95-percent  probability. 

Table  A-1  was  based  on  e  values  near  30  percent,  about  mid  range. 
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TABLE  A-1 


PARAMETER 

REFERRED 

SOURCE 

95.  PERCENT 
CONFIDENCE  LIMIT 

pt'  (P.H.  &  Fuel) 

Data 

±0.6  % 

Pt'  (P.H.  Only) 

Data 

±0.6 

(pt'  P.H.  &  Fuel)  /  (pt'  P.H.  Only) 

Pt '  data 

±1.2 

H, 

Data 

±1.6 

0, 

Data 

±0.8 

^1 

H2  ,  O2 

+  0.7 

y,  +  72 

Ha »  O2  >  yi 

+  7 

V 

Ha,  O2 

+  6 

€ 

J,,  Xa 

+  9 

Tt„ 

Xt 

+  2 

Tt 

+  Xa^’l 

+  5 

Tt/Tto 

Ttg,  Tt 

+  7.5 

Reliability  of  the  Hj  and  0,  analysis  may  be  gaged  by  plotting  fuel^ 
air  equivalence  ratio  for  unburned  fuel  as  calculated  by  Eqs.  (A- 3)  and 
(A-4).  A  plot  of  this  nature  is  shown  in  Fig.  A-1.  An  error  or  deviation 
of  both  measurements  can  cancel  out  or  cause  dispersion  of  data  from  the 
"perfect  agreement"  line. 
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Fig.  A-1  Correlation  of  Fuel  Concentration  from  Hydrogen  and  Oxygen  Analysis 
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APPENDIX  B 

ANALYSIS  OF  CHEMICAL  REACTION  QUENCH  RATES 
IN  THE  GAS  SAMPLING  PROBE 


Since  much  of  the  data  and  conclusions  in  this  report  are  based  on 
sampling  of  the  gas  in  a  burning  stream,  the  question  arises  as  to  just 
how  rapidly  quenching  of  the  chemical  reaction  occurs.  That  quenching 
occurs  is  well  established  by  the  fact  that  combustion  efficiencies,  cal¬ 
culated  based  on  the  gas  analysis,  ranged  from  zero  percent  to  approxi¬ 
mately  90  percent.  The  reaction  is  started  if  any  combustion  can  be 
detected.  The  calculated  combustion  efficiency  indicates  that  the  gas 
entered  the  probe  (Fig.  B-1)  at  some  state  of  reaction  completion,  and 
was  quenched  inside  the  probe  between  the  inlet  orifice  and  a  point  down¬ 
stream,  based  on  the  following  analysis: 

a.  When  the  probe  was  positioned  in  the  supersonic  flow  ahead  of 
the  normal  shock  or  just  aft  of  the  normal  shock  where  no 
emission  was  observed,  the  calculated  combustion  efficiency 
was  zero.  This  was  found  to  be  true  even  though  the  calculated 
combustion  efficiency  one  inch  downstream  was  as  high  as  60 
to  70  percent.  This  observation  indicates  that  quenching 
occurred  faster  than  the  reaction  could  achieve  a  measurable 
start,  or  less  than  3  microseconds  (Fig.  lb). 

b.  When  a  straight  tube  probe  inlet  was  used  instead  of  the  choked 
orifice  of  Fig.  B-  1,  no  quench  occurred  at  all,  and  the  reaction 
proceeded  to  completion  inside  the  probe. 

c.  Analysis  of  the  probe  interior  flow  after  the  gas  had  entered  in¬ 
dicates  that  the  flow  is  choked  and  expands  supersonically.  The 
two  possible  cases  are: 

1.  Separated  flow  (as  shown  in  Fig.  B-1)  where  the  gas  ex¬ 
pands  to  about  M  =  1.  7  or  higher  and  dissipates  into  the 
surrounding  cooler  gas.  (Expansion  and  quenching  occurs 
in  less  than  1  microsecond.  The  flow  is  about  85  percent 
diffused  into  the  surrounding  cool  gases  in  3  microseconds.) 

2.  Fully  developed  expansion,  which  is  possible  with  the  6 

to  7:1  pressure  ratio  available.  (In  this  case,  the  flow  ex¬ 
pands  to  Mach  3.  5,  and  the  reaction  is  quenched  in  less 
than  1  microsecond.  This  flow  shocks  down  in  about 
5  diameters  (Ref.  15).  The  changes  in  static  temperature 
when  passing  the  quenched  gas  through  the  shocks  are  not 
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clear.  Here  again,  if  re-ignition  did  occur,  it  is  expected 
that  its  presence  would  become  known  in  the  gas  analysis. ) 

The  conclusion  reached  from  this  analysis  is  that  quenching  of  the 
reaction  must  occur  within  a  period  of  2  to  3  microseconds. 


30 


STAINLESS  STEEL  GAS  INLET 

0.020  in. 

DIAM— 15  psia 

/ 

— 

-SHOCK  PRESENT 

/ 

y 

IN  SUPERSONIC 

ZONE  ONLY 

\  V  ^ 

^  MIXING  AND 

INR  7nNF 


0.06  in.  DIAM 


COPPER 

WALL  (f('VC'c^ 


WATETmor^ 


Fig.  B-1  Schematic  of  Gas  Sampling  Probe  Quench  Analysis 
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Fig.  1  Supersonic  Combustion  Tunnel  with  Normal  Shock  Configuration 
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Fig.  4  Single  Wedge  Used  to  Produce  Oblique  Shock  Waves 
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Fig.  6  Schematic  of  Pressure  and  Gas  Sampling  Equipment 
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Fig.  7  Mach  Reflected  Normal  Shock  Wave 
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Fig.  8  Variation  of  Flow  Parameters  on  Tunnel  Centerline 
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Fig.  9  Effect  of  Preheater  Temperature  on  Flow 
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Fig.  12  Preheater  Gas  Composition  vs  Total  Temperature 
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Fig.  13  Preheater  Gas  Properties  vs  Total  Temperature 
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b.  Medium  Inlet  Temperature  Level 
Fig,  15  Continued 
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Fig.  16  Combustion  Efficiency  as  a  Function  of  Inlet  Temperature 
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Fig.  17  Variation  of  Flow  Pa 
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Fig.  19  Correlation  of  a  Global  Rate  Equation  for  the  Hydrogen-Oxygen  Reaction 


Fig.  20  Total  Pressure  Loss  vs  Fuel  Concentration  with  Near-Zero  Heat  Release 
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Fig.  21  Total  Pressure  Loss  vs  Fuel  Concentration  at  Various  Combustion  Efficiency  Values 
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Fig.  23  Variation  of  Flow  Parameters  with  Axial  Distance  Aft  of  Normal  Shock  Wave  with  Test  Fuel 
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Fig.  24  Comparison  of  Calculated  Pressure  Loss  with  Experimental  Pressure  Loss  for 
Heat  Addition  Aft  of  the  Normal  Shock 
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Fig.  25  Test  Section  Schlieren  and  Emission  Photograph  with  Test  Fuel  On 
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Fig.  26  Test  Section  Emission  Photograph  with  Test  Fuel  Off 
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